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Abstract 

We characterize all the quasi-stationary distributions and the Q-process associated with a con- 
tinuous state branching process that explodes in finite time. We also provide a rescaling for the 
continuous state branching process conditioned on non-explosion when the branching mechanism is 
regularly varying at 0. 



1 Introduction 

Consider a continuous-state branching process (CSBP for short) t ^ (x) starting from the value x > 
and whose branching mechanism is given by the convex function 

2 i> 

m{^u) = -iu+'^u^ + (e""^ - l + u/il{,,<i})i/(d/i), > (1) 

^ ^(0,oo) 

where 7 G M, a > and 1/ is a Borel measure on (0, 00) such that /^^ 00)^^ ^ h?)v{dh) < 00. It is 
well known that, at any given time i > 0, a; i-> Zt(x) is a subordinator: we denote by u{t, ■) its Laplace 
exponent. In [8], Silverstein showed that the drift dt of this subordinator equals when the CSBP has 
infinite variation paths. The following result gives the expression of this drift when it is not 0. 

Proposition 1.1 When ^ is the Laplace exponent of a Levy process with infinite variation paths, dt = 0. 
Otherwise, for all t > we have 



dt = e-^\ where D := 7 + / h v{dh) (2) 

J(o,i) 
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This entails that dt only depends on the drift of the Levy process whose Laplace exponent is given by ^. 

The goal of this paper is to investigate the quasi-stationary distributions associated with a CSBP that 
explodes in finite time almost surely. From now on, is the distribution on ©([O, +00), [0, +00]) of a 
CSBP starting from x, and Z is the canonical process. More generally, given a probability measure fi on 
[0, 00) we denote by the law of the CSBP initially distributed according to fi. Let T be the lifetime 
of Z, that is, the random time defined by T := Tq A Too where 

(Extinction) Tq := inf{t > : Z^ = 0} 

(Explosion) Too := inf{t > : Z^ = 00} 

By quasi-stationary distribution (QSD for short), we mean a probability measure /x on (0, 00) such that 

F^{Zte-\T >t) = 

In that case, it is a simple matter to check that under F^ the random variable T has an exponential 
distribution, the parameter of which is called the rate of decay of /x. 

A brief review of the litterature: the extinction case. 

Li [7] and Lambert [^] considered the extinction case T = Tq < 00 almost surely and studied the 
CSBP conditioned on non-extinction. Let us recall some of their results. When ^ is subcritical, that is 
\1''(0) > 0, there exists a family (/i^; < /? < ^''(0)) of quasi-stationary distributions on (0, 00) where 
/3 is the rate of decay of j^Lf,. These distributions ai^e characterized by their Laplace transfomis as follows 

/ n^idr)e-''^ = 1 - e-'^*^^), VA > (3) 

J(0,oo) 

where $(A) := /^°° du/^{u). For any (3 > ^'(0) they proved that there is no QSD with rate of 
decay /3, and that Equation (3) does not define the Laplace transform of a probability measure on (0, 00). 
Additionally, the value /3 = 'I''(0) yields the so-called Yaglom limit: 

F^{Zt e-\T>t) ^ ^%'(o)(-), Vx>0 

When ^ is critical, that is ^''(0) = 0, the preceding quantity converges to a trivial limit for all 2; > and 
Equation (3) does not define the Laplace transform of a probability measure on (0, 00). However, under 
the condition \1'"(0) < 00, they proved the following convergence 

F.(|>.|T>i)^-.exp(-^), Vx>0,.>0 (4) 

In Section 3, we will provide a generahsation of this result only assuming that ^ is regularly varying at 
0, see the forthcoming Proposition 3.1. 

Finally in both critical and subcritical cases, for any given value t > the process {Zr,r G [0, t]) 
conditioned on s < T admits a limiting distribution as s — > 00, called Q-process. The law of the 
Q-process is obtained as a /i-transform of P as follows 

dQx\Tt ■■= d^x\Tt, Vx > 
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Main results: the explosive case. 

Suppose T = Too < oo almost surely and remark that conditioning a CSBP on non-explosion does not 
affect the branching property. Subsequently the law of conditioned on T > f is infinitely divisible: if 
it admits a hmit as t goes to oo, the limit has to be infinitely divisible as well. 
In [ ], Grey proved that the almost sure explosion is equivalent with 

< 0, Vu > and / — ^ < oo (5) 

In particular, the first condition ensures that is the Laplace exponent of a subordinator. We introduce 
the notation: 

'I'(+oo) := lim ^{u) G [-oo, 0) 

u~^oo 

It is a simple matter to verify that 

^'(+oo) = -oo u{0, oo) = oo or < 

Conversely, ^(+oo) G (— oo, 0) imphes ^(+oo) = — i^(0,oo). Our result below shows that this quantity 
plays a role analogue to ^''(0+) in the extinction case. 

Theorem 1 Suppose T = T^o < oo almost surely and set 

For any f3 > there exists a unique quasi- stationary distribution /i^ associated to the rate of decay /3 
which is infinitely divisible and is characterized by 

f M^(dr)e-"^ = e-'^*W, VA>0 (6) 

Additionally, the following dichotomy holds true: 

(i) ^(+oo) G (—00,0). The limiting conditional distribution is given by 

lim P^(Zt G • I T > t) = ^f.,(o.^)(-), Vx G (0, oo) 

(ii) ^'(+oo) = — oo. The limiting conditional distribution is trivial: 

lim P^(Zf < a 1 T > t) = 0, Va, X G (0, oo) 

Let us make some comments. Firstly this theorem implies that A i-^- $(A) is the Laplace exponent of 
a subordinator, and so, is the distribution of a <I>-Levy process taken at time /3. Secondly there is a 
similarity with the extinction case: the limiting conditional distribution is trivial i.f.f. ^(+oo) = — oo 
so that the dichotomy on the value ^'(+oo) is the explosive counterpart of the dichotomy on the value 
^'(0) in the extinction case. Also, note the similarity in the definition of the Laplace transforms (3) and 
(6). However, there are two major differences with the extinction case: firstly there is no restriction on 
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the rates of decay. Secondly, even if the limiting conditional distribution is trivial when ^'(+oo) = — oo, 
there exists a family of QSDs. We now propose, when is regularly varying at 0, a rescaling of the 
CSBP conditioned on non-explosion that ensures a non-trivial limit. Recall that we call slowly varying 
function at any continuous map L : (0, oo) — > (0, oo) such that for any a S (0, oo), L{au) / L{u) — > 1 
as u I 0. 

Theorem 2 Suppose that ^{u) = —v}-~"^L{u) with L a slowly varying function at and a G (0, 1), and 
assume that ^'(+oo) = — oo. Consider any function f : [0, oo) — )• (0, oo) satisfying ^(/(t)~^)/(t) ~ 
^{u{t, 0+)) as t ^ oo. Then the following convergence holds true: 



^xzt/fit) U < T 



t^oo 



^-xXya^ Vx,Ag(0,Oo) 



Notice that such a function / exists. Indeed, the map u i-^ ^{u)/u is strictly increasing and takes all 
values in (— oo, D), the existence then follows from simple arguments. 

Example 1.2 When ^'(n) = -ku^''^ with k > and a £ (0, 1), we have f{t) ~ (a/ct)^!-")/"'' as 
t oo. When ^(u) = -cu-ku^-°' with k,c> Oanda £ (0, 1), we have f{t) ~ (A;/c)(^-")/"'e'=*/" 
as t ^ oo. 

Finally, we compute the Q-process associated with an explosive CSBP. Let Tt be the sigma-field gener- 
ated by {Zs, s G [0, t]), for any t G [0, oo). 

Theorem 3 We only assume that T = Too < oo almost surely. For each x > 0, there exists a distribu- 
tion on ID)([0, oo), [0, oo)) such that for any t > and any event G J-f 

lim F,(e|T > s) = Q,(e) 

S— >-CXD 

Furthermore, Q_x is the law of the ^^-CSBP where 

^!^{u) = Du 

The Q-process appears as the ^'-CSBP from which one has removed all the jumps: only the deterministic 
part remains (recall the expression of the drift from Proposition 1.1). Notice that the Q-process cannot be 
defined through a /i-transform of the CSBP: actually the distribution of the Q-process on B([0, t], [0, oo)) 
is not even absolutely continuous with respect to that of the ^-CSBP, except when the Levy measure u 
is finite. 



The paper is organized as follows. In the second section, we examine the discrete case: as no 
explosion occurs in discrete time, we need to consider continuous-time Galton- Watson processes. We 
provide a complete description of the QSDs when this process explodes in finite time almost surely. In 
the third section, we recall basic facts about CSBPs then we prove Propositions 1.1 and 3.1. In the fourth 
section, we prove Theorems 1, 2 and 3. 
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2 The discrete case 



A discrete-state branching process (2^, t > 0) is a continuous-time Markov process taking values in Z+ 
that verifies the branching property (see for instance [ 1 ], Chapter III). It can be seen as a Galton-Watson 
process with offspring distribution ^ where each individual has an independent exponential lifetime with 
parameter c > 0. Let us denote by (/)(A) = '^'^Lo ^^^{^)^ ^ [Oil] the generating function of 
the Galton-Watson process. We denote by P„ the law on the space ©([O, oo), Z+) of Z starting from 
n G Z+, and E.„ the related expectation operator. More generally, we denote by the law of Z 
with initial distribution fi, where /i is a probability distribution on Z+. The semigroup of the DSBP is 
characterized via the Laplace transform 

VrG (0,l),Vt e [0,oo), E„ [ r^* ] =F(i,rr where j = t (7) 

We denote by t the lifetime of Z, that is, the infimum of the extinction time Tq and the explosion time 
Too: it is well-known that P„(to < t) = F{t, 0+)" and P„(Too < t) = 1 - F{t, 1-)". In this section, 
we assume that there is explosion in finite time almost surely: this implies that the smallest solution to 
the equation = x equals and that Jj^_ < oo. We define 

r dx 
J I c((/)(x)-x) 

Cleaiiy r i-> <I>(r) is the inverse function of t ^ F{t, 1-), that is for all t > 0,^{F{t, 1-)) = t. We 
say that a measure /i on N = {1, 2, . . .} is a quasi-stationary distribution (QSD) for Z if 

P^{Zt G ■\r>t) = fi{-) 

From the Markov property, we deduce that t has an exponential distribution under P^, the parameter of 
which is called the rate of decay of fi. 

Theorem 4 Suppose there is explosion infinite time almost surely. Let (3q := c(l — There is a 

unique quasi-stationary distribution /i^ associated with the rate of decay /3 if and only if P is of the form 
n/3o, with n € N. /f is characterized by its Laplace transform 

J]M,s(A^)^' = e-''*(^\ VrG(0,l] 

k 

For any initial condition n € N we have 

lim P.„( G • I T > t ) = /^n8o(-) 

I— >oo 

Finally the Q-process is the constant process. More precisely denote by Ft the sigma-field generated by 
{Zg, s G [0, t])for any t > 0. For any event Q G Ft and any initial condition n G N 

lim p„(e|T>s) = Q„(e) 

where Q is the distribution of the DSBP with the trivial generating function F{t, r) = r. 
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Let us make some comments. First there exists only a countable family of QSDs. This is due to the 
restrictive condition that our process takes values in N. Also, observe the similarity with Theorems 1 
and 3: indeed a DSBP can be seen as a particular CSBP starting from an integer and whose branching 
mechanism is the Laplace transform of a compound Poisson process with integer-valued jumps. In 
particular = c^(/c + 1) for every integer k > 1. Hence the quantity c(l — ^(1)) in the DSBP 

case corresponds to v{0, oo) in the CSBP case. 

Sketch of the proof. We only prove that /3 has to be of the fomi Pqu, with n G N. The other properties 
derive from simple calculations. Fix n G N and observe that 

E„[r^* \ i>t\ — )■ exp(-^>(r)n/3o) 
From this last convergence and the fact that <I>(r) — > 0, we deduce that r i— )■ exp(— $(r)n/3o) is the 

rtl 

Laplace transform of a probability measure on Z+. As ^{r) — > +oo, we deduce that this probability 

measure does not charge 0. Also, a simple calculation ensures that //^^({l}) > 0. 

Conversely, suppose that r i-> e~^^^'^^ is the Laplace transform of a probability measure on N. Denote 

by m G N the smallest integer such that > 0. Then we have for all r G (0, 1) 



/^;3o({l}> + E/^^o(W) 



A;>1 



This implies that 



and so, m = G 



/x,({m})r- ~ (^,„({l})r)/o 



3 Continuous-state branching processes, drift and critical case 

Recall the form taken by a branching mechanism ^ from (1). A continuous-state branching process 
associated with is a Markov process whose semigroup is characterized by the following identity 

E^[e~^^'] = e-^"(*'^\ VA > (8) 

where 

dtuit, A) = -^[u{t, A)), n(0, A) = A (9) 

We will denote by q the second root of ^ that is q := sup{n > : ^(u) < 0}. Note that g is in the 
subcritical and critical cases {i.e. when ^''(0) > 0) while it equals oo when ^ is the Laplace exponent of 
a subordinator. 

Classical results [3] ensure that A i— u{t, A) is the Laplace exponent of a subordinator. Thanks to the 
Levy-Khintchine formula, there exist at,dt > and a Borel measure Wf on (0, oo) with /(Qf^)(l A 
h)wt{dh) < oo such that for all A G (0, oo) 

u{t,X) = at + dtX+ [ {l-e-^^)wtidh) (10) 

^(0,oo) 
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3.1 Proof of Proposition 1.1 

For dt to be positive it is necessary and sufficient that 

fj = and / {1 Ah) v{dh) < oo 

i{0,oo) 

see Corollary in [8]. This is equivalent to saying that ^' is the Laplace exponent of a Levy process with 
finite variation paths. We suppose this condition fulfilled until the end of the proof. In this case, we can 
write 



^{u) = Du+ [ {e-""^ - 1) u{dh) 

J(O.oo) 



'{0,oo) 

Results [2] on the Laplace exponent of a Levy process with finite variation paths ensure that 

— > ^ , — ^ — ^ — > dt 

Also, from the very definition of u{t, A) one immediately gets for all t > 0, A > 

^ du 



(11) 



lu(t,X) ^(-U) 

Now remark that u{t, A) ^ oo as A ^ oo and write for any t > 0, A > 

\ X J Jq u{s,X) Jq u{s,X) 

We have ^{u)/u — > D as u — )■ oo and s ^ u{s, A) is monotone: the dominated convergence theorem 
implies that log{u{t, A)/A) — )• —Dt as A — > oo. This ends the proof. ■ 

3.2 A more general result in the critical case with extinction 

In [9], Slack showed that any critical Galton- Watson process with a regularly varying generating function 
can be properly rescaled so that, conditioned on non-extinction, it converges towards a non-trivial limit. 
In the present paper, we prove a similar result for critical CSBPs that become extinct almost surely. 

Proposition 3.1 Suppose that ^(n) = ti^~'""L(u) with L a slowly varymg function at and a S (0, 1]. 
Assume that T = Tq < oo almost surely. Fix any function f : [0, oo) (0, oo) verifying f{t) ^ u{t, oo) 
as t ^ oo. Then we have the following convergence 

Eje-^^'^W|t<Tl ^ 1- (i + A-")-'/", Vx,AG(0,oo) 

We recover in particular the finite variance case (4) of Lambert and Li. Our result also covers the so- 
called stable branching mechanisms ^'(ti) = u^"*^" with a G (0, 1]. Once again one should notice the 
similarity between the explosion and the extinction cases: the assumption made in this proposition on 
the branching mechanism is quite similar to that of Theorem 2. 
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Proof The proof is similar to that of Theorem 1 in [9], so that some calculations are left to the reader. 
Recall that the extinction in finite time implies that for any t > 0, u{t, oo) < oo. Fix x, A € (0, oo). 
First, remark that for any r E (0, oo) we have 
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1 



r u(r, oo)"L(ii(r, oo)) r Jq 



^'(n(s,oo)) - n(s,oo)^'(u(s,oo)) 
OsU[s, OO) : — as 



^'(n(s,oo))' 

Since ^ is regularly varying at we use Theorem 2 in [6] to deduce that tt*' (n) (n) — )■ 1 + a as u I 0. 
Hence we get 



u(r, oo)"L(n(r, oo)) ~ as r oo 

a r 



(12) 



Since r i-^ u{r, oo) decreases from oo to 0, for any t E (0, oo) there exists s{t) = s G (0, oo) such that 
u{s, oo) = Xf{t). Since L is slowly varying and f{t) ~ u{t, oo) as i — )• oo, we use (12) to deduce that 
A°s ~ t as t — )• oo. Also notice that 



log 



u{t + s, oo) 
u(t, oo) 



dru{r, oo) 
u{r, oo) 



dr 



t+s 



u{r, oo)" L[u{r, oo))dr 



~ log(l + A 

t->oo a 



where we use (12) at second line. Then, we have 



~ 1 

t—^oo 



^-xu{t,\f{t)) _ ^-xu{t,oo) 
g— XM(i,oo) 

u(t + s, oo) 



u{t, oo) 
1-(1 + A-")-^/" 



This ends the proof. 



4 Quasi-stationary distributions and Q-process in the explosive case 

In this section, we consider a branching mechanism ^ which verifies (5) so that Z explodes in finite time 
almost surely. An elementary calculation entails for all i > and all x > 

P^(T >t) = e-^"(*'0+) 

We introduce 

This increasing function admits a continuous inverse, namely the function t i-)- u{t, 0+). Also, thanks to 
Equation (11) we deduce the identities 

<^{u{t,X)) = t + <l>{X) , u{t,X) = u{t + <^{X),0+) , Vi,A>0 (13) 

Finally we introduce the quantity ^(+oo) := lim ^(u) G [— oo,0). The dominated convergence 
theorem ensures that ^'(+oo) = — z^(0, oo) as soon as ^(+oo) is not infinite. 
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4.1 Proof of Theorem 1 

First we compute the necessary form of the QSD. Fix /3 > 0. We get for all t > 

e-/5* = r„jT >t)= f M«((ir)e-^"(*'0+) 



P^,(T>t)= / ^i,{d 

J(O.oo) 

Letting t = <&(A) for any A > we obtain 



i(0,oo) 

Consequently there is at most one QSD corresponding to the rate of decay (3. Now suppose that /i^ is a 
probability distribution on (0, oo) then the following calculation ensures that it is quasi-stationary. For 

any A € (0, oo) we have 



P^,(T>t) PM,,(T>t) 



/(o,oo) (f^^')e-" "(*'°+) e-/3 *('^(*.o+)) 

We now assume 'I'(+oo) e (— oo, 0) and we prove that A i— > e~^*^^^ is indeed the Laplace transform of 
a probability measure /i^ on (0, oo). Let x := P/v{0, oo), for all A > we have 

[e-^^* I T > t] = "^jp^^^'^^^ ^ = exp ( - x{u{t + cI>(A), 0+) - u{t, 0+)) 
From Equation (1 1) we get that u{t, 0+) — ^ oo as t — ^ oo and that 

A{t+${A),0+) ^l^J 

Since ^'(n) — )• — j^(0, oo) as li — )■ oo, one deduces that 

^"(*'°+) (iu n(t + $(A),0+) -n(t,0+) 



„(j+$(A),0+) ^(^) 1^(0, oo) 



and therefore 

-AZt|rp ^ ^ p-$(A)xi/(0,oo) 

4— >oo 



E^. [e i > t| — > e 



Since <&(A) ^> as A J, 0, we deduce that A i-)' g-*{A) 2:!y{o,oo) j-j^g Laplace transform of a probability 
measure on [0, oo). Moreover, it does not charge since <I>(A)— )-ooasA^oo. 

We now suppose ^(— oo) = — oo. An easy adaptation of the preceding arguments ensure that for any 

X, A > 

Ea;[e~^^*|T >t] — ^ 
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Hence the limiting distribution is trivial: it is a Dirac mass at infinity. However, let us prove that fis is a 
well-defined probability distribution on (0, oo). For every e > 0, define the branching mechanism 

^e{u) := [ {e-^'' -l){iy{dh) + -6^ne{dh)) 

J(e,oo) e 



-(e^^"-l)+ / (e"'^" - 1) 
e J(t.oo) 



'(e,oo) 

Observe that for any n > 0, ^'e(n) J, ^'(n) as e J, 0. Thus by monotone convergence we deduce that 

The first part of the proof applies to and therefore the l.h.s. of the preceding equation is the Laplace 
transform taken at A of an infinitely divisible distribution on (0, oo). Since the limit is continuous at 
and vanishes at oo, we deduce that it is the Laplace transform of an infinitely divisible distribution on 
(0, oo). The asserted property follows. ■ 

4.2 Proof of Theorem! 

Fix A, X € (0, oo). For any t € (0, oo), we have 

-ilogE,.[e-^^'/^W |t < T] = u{t,Xf{t)~^)-u{t,0+) 
The proof relies on two lemmas, whose proofs ai^e postponed to the end of the subsection. 
Lemma 4.1 As u 10, we have <^(n) ~ u/{—a^{u)). 
Observe that f{t) — > +oo as t ^ oo. The lemma now implies 

Since L is slowly varying at 0+, we deduce that ^'(A/(t)^^) ^ A^~"^'(/(t)~^) as t — > oo. Thus the 
very definition of / entails 

*(M(t,0+))$(A/(t)-^) ~ -A"a-^ (14) 



Lemma 4.2 The following holds true as t ^ oo 

ru(tfi+) fu{t,0+) 



lu{t,x /{<)-!) ^(v) Ju{t,x /{t)-i) ^iu{t, 0+)) 
From the latter lemma, we deduce 



u{t,Xf{t)-^)-u{t,0+) ~ -^{u{t,0+)) 



«(t,A/(i)-i) 



~ -M/(n(t,0+))<&(A/(t)-i) 
^ A a 
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where we use Equations (11) and (14) at the second and third line. The theorem is proved. 

Proof of Lemma 4.1. Recall the definition of An integration by parts yields that for all u E [0, oo) 

$(n) = — + / — r^dv 



Recall from Theorem 2 in [6] that v"^' {v) /"^{v) —t- 1 — a as v J, 0. Therefore an elementary calculation 
ends the proof. ■ 
Proof of Lemma 4.2. For alH G [0, cx)) we have 

A(i,A/(i)-i) ~ Ju{t,\f(t)-^) "^RW)) ~ Ju(t,xm-') ~Mv)Mu{t^^W' 
The convexity of implies that for all v € [u{t, 0+), u(t, A f{t)~^)] we have 

< ^(n(t, 0+)) - ^(v) < -^'{u{t, 0+)) [v - u{t, 0+)) (15) 

Since '^'{v) /"^{v) goes to as u ^> oo, the following identity (and therefore the lemma) 

u{t,o+) ,^fu(t,0+))-^(v) , . rHt,o+) dv 

dv — ' 



would be ensured by (15) and the boundedness of t u{t, A f{t)~^) — u{t, 0+), that we now prove. 
Fix k G {—2D, oo). Since "^'{v) converges to D as w — > oo, for t large enough we get from (15) that 
^'(ii) > 'i!{u{t,0+))-k{v-u{t,0+))foTallv £ [u{t,0+),u{t, A /(t)"^)]. A simple calculation then 
yields 

Using log(l + v) > v/2 for v small, and since <l>(A/(t)~'^) — )• 0, we get for t large enough 

2^('u(t,0+)) - ^ ■'^ ^ ' 

From Equation (14), we deduce that t ^ u{t, A f{t)~^) — u{t, 0+) is bounded. ■ 

4.3 Proof of Theorem 3 

Fix t > 0. Since we are dealing with non-decreasing processes and since the asserted limiting process 
is continuous, the convergence of finite-dimensional marginals suffices to prove the theorem. Hence we 
will prove that for any n > 1, any n-uplets 0<ti<...<tri<i and any coefficients Ai, . . . , A^ > 
we have 

lim --log(E^[e-^i^*i---^"^*"|T>t + s]) = Aidt, + ... + A„dt„ (16) 

s— >oo X 

Thanks to an easy recursion, we get 

-- log(E^[e~^i^'i~-^^"^'" IT > t + si) 

X 

= U^i, Ai + u(^t2 - ti, X2 + ■ ■ ■ + u(tn - tn-1, Xn + u{t + S - tn,0+)) • • • )^ - u{t + S, 0+) 



11 



To prove (16), we proceed via a recun^ence on n. We check the case n = 1. We know that A — )■ u{t, A) 
is a concave function from to IR+ such that 

— ^ as A — ;> oo 

A 

The concavity of A — )■ u(t, A) then implies that d\u{t, A) — ?> as A — )■ oo. Writing u{t + s, 0+) = 
u{ti,u{t + s — ti, 0+)) one deduces that 

u[ti, Ai + u{t + s — ti, 0+)) — u{t + s, 0+) —7- Aifij^ as s — ;> oo 

The result forn = 1 follows. The end of the recursion is elementaiy. ■ 
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